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1.  Background 

The  basic  intent  of  this  research  was  three-fold: 

(1)  to  implement  a  soil  hydrology  model  (BATS)  in  a  mesoscale  atmospheric  model, 
specifically  the  Penn  State  Mesoscale  Model,  version  5  (MM5); 

(2)  to  determine  the  best  way  to  arrive  at  a  set  of  optimum  initial  soil  water  content 
values  using  a  land  surface  process  scheme  implanted  in  MM5  called  the  Biosphere- 
Atmosphere-Transfer-Scheme  (BATS); 

(3)  to  study  a  phenomenon  we  call  ‘decoupling’  or  rapid  soil  drying  (Capehart  and 
Carlson,  1997).  Our  original  idea  was  to  invert  the  results  of  the  MM5/BATS  package 
in  order  to  determine  optimum  values  of  soil  water  content.  The  latter  proved 
untenable. 

Optimum  soil  water  content  values  in  a  land  surface  process  model  can  be  defined  as  thos 
value  of  soil  moisture  which  will  result  in  a  minimization  of  errors  in  afternoon 
temperature  and  humidity  forecasts  using  MM5.  By  rapid  soil  drying  we  refer  to  a 
phenomenon  discussed  by  Capehart  and  Carlson  (1997),  where  by  the  surface  of  the  sc  l 
dries  out  very  rapidly  (in  a  few  hours)  leaving  a  relatively  wet  sub  soil.  Rapid  soil  dr;,;iig 
is  accompanied  by  a  rise  in  surface  air  temperature  and  a  decrease  in  relative  and  aboluste 
humidity  due  to  a  decrease  in  evapotranspiration  and  an  increase  in  sensible  heat  flux 
from  the  ground. 

Initially,  the  land  surface  component  in  MM5  called  the  Biosphere/ Atmosphere/Transfer 
Model  proved  inadequate  to  study  rapid  soil  drying  for  two  reasons.  First,  the  originally 
proposed  method  for  obtaining  optimum  surface  parameters  proved  unwieldy  and 
unreliable.  Second,  BATS  lacked  sufficient  vertical  resolution.  Our  research  has  shown 
that  rapid  soil  drying  occurs  in  just  the  two  1  or  2  cm  of  the  soil.  Predicting  the  onset  of 
rapid  soil  drying  is  dependent  on  knowledge  of  both  a  thin  surface  layer  as  well  as  the 


soil  water  content  in  deeper  layers.  It  is  also  highly  dependent  on  soil  type  and 
vegetation,  fraction. 

Thus,  in  the  middle  of  the  granting  period  it  was  necessary  for  us  to  rethink  the  problem 
from  scratch  and  to  begin  again.  Fortunately,  a  more  highly  resolved  version  of  BATS 
(called  SHEELS)  was  being  developed  at  that  time  by  scientists  working  at  NASA, 
Marshall  Space  Flight  Center,  Huntsville,  AL.  Accordingly,  we  used  SHEELS  as  the  land 
surface  component  in  MM5  and  a  soil  hydrology  model  (called  SHM)  to  initialize  the  soil 
water  content  was  developed  during  an  earlier  phase  of  our  contractual  relationship  with 
the  Air  Force. 

Fractional  vegetation  cover  proved  also  to  be  of  great  importance  in  predicting  afternoon 
temperature  and  humidity.  Realizing  the  importance  of  fractional  vegetation  cover,  we 
created  a  data  base  on  the  web  of  fractional  vegetation  cover  data  for  seven  years  (1990- 
1997)  at  1  km  resolution  for  part  of  North  America  from  AVHRR  satellite  data;  this  data 
is  now  available  on  our  web  site  referred  to  later  in  this  report.  Ultimately  we  were  able 
to  take  these  measured  values  of  fractional  vegetation  cover  derived  from  satellite 
imagery  and  soil  moisture  values  determined  with  a  soil  hydrology  model  to  initialize  the 
MM4/SHEELS  package.  We  then  made  a  series  of  simulations  with  MM5/SHEELS  for 
one  case  study  and  subsequently  verified  the  output  of  these  simulations  against  surface 
weather  and  soil  data  taken  in  at  selected  points  in  the  Oklahoma  meso  network. 

For  details  of  this  research  project,  the  reader  is  referred  to  the  attached  paper  by 
Santanello  and  Carlson  (2001);  copies  are  attached  to  this  report.  (A  parallel  paper 
published  by  scientists  at  NOAA,  NSSL  in  Norman  OK  (Crawford  et  al.,  2000)  was 
inspired  by  this  research.) 

2.  Procedures 

Two  important  parameters  are  required  to  specify  soil  moisture  behavior  in  land  surface 
components  such  as  SHEELS:  fractional  vegetation  cover  and  soil  moisture  content.  The 


vegetation  fraction,  with  a  default  soil  water  content  and  a  calculated  (from  NDVI) 
vegetation  fraction,  and  with  a  calculated  soil  water  content  and  a  calculated  vegetation 
fraction.  Simulations  were  performed  for  the  middle  third  of  the  United  States,  but 
centered  over  the  Oklahoma  meso-network.  The  domain  consisted  of  three  nested 
domains,  a  coarse  outer  grid  of  36-km  spacing,  an  interior  grid  of  12  km,  and  a  fine  mesh 
of  4  km.  The  Utter  corresponded  to  the  meso-network  itself. 

The  simulations  were  chosen  for  16  July,  1996,  a  period  of  abundant  sunshine  and  light 
winds  immediately  following  a  moderate  precipitation  event.  We  chose  these  conditions 
because  of  the  likelihood  that  rapid  soil  drying  would  occur.  Simulations  were  made  out 
to  36  hours  starting  at  1200  UTM,  in  order  to  show  the  daytime  rise  in  temperature  on 
two  successive  days.  Simulations  had  to  be  run  several  times  before  it  was  felt  that  the 
model  was  behaving  as  it  should.  Following  each  of  the  first  two  or  three  series  of 
simulations  we  discovered  errors  in  SHEELS  which  required  an  extended  period  of 
remediation  lasting  almost  a  year  during  which  time  the  creators  of  SHEELS  at  NASA, 
Huntsville,  were  able  to  produce  successively  improved  versions  of  the  model. 
Altogether,  the  duration  required  to  determine  that  BATS  was  inadequate  for  the  job  and 
the  period  of  remediation  at  NASA  for  SHEELS  lasted  well  over  one  year. 

3.  Results 

Results  showed  that  rapid  soil  drying  can  profoundly  affect  afternoon  temperature  and 
humidity,  but  it’s  spatial  distribution  is  necessarily  quite  variable.  However,  rapid  soil 
drying  depends  on  both  soil  water  content  and  fractional  vegetation  cover  as  well  as  on 
soil  type  and  cloud  cover.  We  found  that  rapid  soil  drying,  with  a  concomitant  rise  in 
sensible  heat  flux  and  a  decrease  in  evapotranspiration,  depended  on  a  rather  critical 
threshold  of  soil  water  content,  soil  type  and  vegetation  cover.  Rapid  soil  drying  would 
occur,  moreover,  only  within  a  relatively  narrow  range  of  values.  For  example,  rapid  soil 
drying  would  occur  when  the  average  soil  water  content  over  the  top  10  cm  of  the  soil 
began  to  fall  below  about  0.6  of  field  capacity;  this  rapid  drying  would  occur  only  until 
the  soil  water  content  reached  about  0.3  of  field  capacity.  It  would  occur,  moreover,  only 


when  the  fractional  vegetation  cover  was  below  about  0.6.  These  threshold  values  varied 
considerably  with  soil  type  (soil  hydraulic  conductivity).  For  that  reason  the  simulations 
showed  a  highly  variable  distribution  of  surface  sensible  heat  flux  and  evapotranspiration 
in  response  to  the  variations  in  the  various  land  surface  properties.  A  table  of  threshold 
values  is  presented  in  the  paper  by  Santanello  and  Carlson  (Table  3). 

Because  of  a  problem  in  translating  the  soil  hydrology  model  profiles  to  SHEELS  soil 
water  content  profiles,  the  hydrology  model  somewhat  overestimated  the  soil  water 
content  for  the  initial  simulation  time  of  the  MM5/SHEELS.  As  a  result,  the  initial  soil 
water  content  from  the  hydrology  model  was  above  the  threshold  value  for  rapid  soil 
drying.  Simulations  with  the  valued  derived  from  the  hydrology  model  thus  produced 
rather  incorrect  afternoon  temperatures  that  were  too  low  when  validated  against  ground 
observations.  As  the  default  soil  water  content  values  were  below  the  threshold, 
simulations  with  the  default  soil  water  content  values  proved  more  accurate,  as  these 
simulations  were  able  to  capture  the  rapid  soil  drying.  However,  it  is  our  feeling  that, 
once  the  problem  of  converting  the  soil  hydrology  model’s  profile  to  that  of  SHEELS  is 
solved,  more  accurate  patterns  of  soil  water  content  would  be  produced  by  the  SHM  than 
for  the  default  case.  The  allied  paper  by  Crawford  et  al.,  (2000)  also  makes  the  same 
suggestion. 

The  vegetation  cover  analyses  pro’-^d  to  be  more  useful  than  the  soil  water  content,  as 
rapid  soil  drying  is  very  sensitive  to  the  amount  of  vegetation  cover.  Simulations  of 
afternoon  temperature  and  humidity  were  much  less  accurate  when  the  default  vegetation 
fields  were  used  simply  because  the  latter  were  much  more  uniform  than  the  actual 
pattern  of  vegetation. 

4.  Conclusions  and  Recommendations 


We  were  able  to  arrive  at  a  few  important  conclusions. 


•  Fractional  vegetation  cover  derived  from  biweekly  remotely  sensed  NDVI  data 
provides  some  improvement  over  default  estimates. 

•  Soil  moisture  simulated  from  a  hydrological  model  with  finely  resolved  vertical 
resolution  shows  signs  of  rapid  soil  drying  at  the  surface  under  appropriate 
atmosphere  and  soil  conditions. 

•  Higher  vertical  resolution  soil  models,  such  as  that  used  ’’i  SHEELS,  have  the  ability 
to  simulation  rapid  soil  drying.  However,  coarser  models  do  not.  We  recommend  that 
all  land  surface  process  models  resolve  the  top  1  cm  of  soil  layer  in  order  to 
accurately  predict  rapid  soil  drying  and  its  consequences. 

•  Rapid  soil  drying  occurs  within  the  top  1-2  cm  layer  when  the  average  soil  water 
content  over  the  top  10  cm  falls  into  a  certain  range  of  values  that  are  dependent  on 
soil  type.  It  is  not  necessary  to  accurately  specify  the  soil  water  content  in  the  top  1  - 
2  cm  because  a  good  land  surface  scheme  will  quickly  reproduce  the  correct  drying. 

•  The  effect  of  near-surface  soil  water  content  and  soil  type  is  much  larger  than  that  of 
Fr  on  soil  water  content  in  the  surface  layer,  with  the  exception  of  regions  where  the 
latter  exceeds  0.8.  More  accurate  specification  of  fractional  vegetation  cover  is 
required  for  more  accurate  predictions  of  afternoon  temperature  and  humidity. 

•  More  realistic  soil  water  content  values,  as  derived  from  a  hydrological  model,  might 
improve  afternoon  temperature  and  humidity  forecasts,  although  this  was  not 
demonstrated  in  this  study. 
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ABSTRACT 

Rapid  soil-surface  drying,  which  is  called  “decoupling,”  accompanied  by  an  increase  in  near-surface  air 
temperature  and  sensible  heat  flux,  is  typically  confined  to  the  top  1-2  cm  of  the  soil,  while  the  deeper  layers 
remain  relatively  moist.  Because  decoupling  depends  also  on  a  precise  knowledge  of  fractional  vegetation  cover, 
soil  properties,  and  soil  water  content,  an  accurate  knowledge  of  these  parameters  is  essential  for  making  good 
predictions  of  temperature  and  humidity.  Accordingly,  some  simulations  centered  on  the  Atmospheric  Radiation 
Measurement  Program  Cloud  and  Radiation  Test  Bed  Southern  Great  Plains  site  in  Kansas  and  Oklahoma  using 
a  high-resolution  substrate  layer  (Simulator  for  Hydrology  and  Energy  Exchange  at  the  Land  Surface),  the  Fifth- 
Generation  Pennsylvania  State  University-National  Center  for  Atmospheric  Research  Mesoscale  Model,  and 
derived  and  default  values  for  soil  water  content  and  fractional  vegetation  cover  are  presented.  In  so  doing,  the 
following  points  are  made:  1)  decoupling  occurs  only  within  certain  threshold  ranges  of  soil  water  content  that 
are  closely  related  to  the  soil  type  and  2)  a  knowledge  of  fractional  vegetation  cover  derived  from  concurrent 
observations  is  necessary  for  capturing  the  spatial  variation  in  rapid  soil  drying  in  forecast  models. 


1.  Introduction 

The  interaction  between  the  land  surface  and  atmo¬ 
sphere  is  crucial  for  predicting  daily  maximum  air  tem¬ 
perature  and  humidity.  Errors  in  predicting  these  vari¬ 
ables  result,  in  part,  from  incorrectly  determining  the 
partition  of  solar  radiation  into  surface  turbulent  fluxes 
of  heat  and  moisture.  Studies  have  shown  that  the  two 
most  important  land  surface  parameters  that  influence 
the  surface  energy  budget  are  the  near-surface  soil  water 
content  (SWC)  and  the  fraction  of  vegetation  cover  (Fr ; 
Gillies  et  al.  1997;  Sun  and  Bosilovich  1996).  SWC  is 
the  most  important  soil  characteristic  governing  fluxes 
at  the  air-soil  interface,  its  impacts  being  much  greater 
than  those  of  soil  albedo  and  soil  type  (Smith  et  al. 
1994;  McCumber  and  Pielke  1981).  Moreover,  surface 
heat  fluxes  are  one  order  of  magnitude  more  sensitive 
to  the  SWC  profile  than  they  are  to  the  soil  temperature 
profile  (McCumber  and  Pielke  1981).  SWC  also  plays 
a  role  in  determining  the  heat  and  moisture  fluxes  within 
the  soil  (McCumber  and  Pielke  1981).  The  ability  of 
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new  land  surface  schemes  to  handle  more  complex  rep¬ 
resentations  of  the  land  surface  highlights  the  need  for 
more  accurate  estimates  of  SWC  (and  its  vertical  dis¬ 
tribution). 

Despite  its  importance  in  atmospheric  forecast  mod¬ 
els,  SWC  still  remains  an  elusive  variable  for  real-time 
weather  prediction  because  of  the  lack  of  field  mea¬ 
surements,  limited  accuracy  of  remote  sensing  methods, 
and  the  overall  heterogeneous  nature  of  soil  properties 
that  cannot  be  captured  by  models  with  grid  spacing  on 
the  order  of  kilometers  (Capehart  and  Carlson  1997, 
hereinafter  CC97;  Capehart  1996).  Soil  moisture  esti¬ 
mates  made  from  space  using  microwave  and  thermal 
radiometry  show  some  skill  in  assessing  surface  mois¬ 
ture  on  large  spatial  scales,  but  these  techniques  fail  to 
capture  the  vertical  variability  in  rapid  soil  drying  that 
can  occur  as  a  result  of  a  process  discussed  by  CC97 
called  “decoupling.” 

Decoupling  occurs  when  rapid  soil  drying  takes  place 
under  intense  sunshine,  leading  to  a  marked  reduction 
in  soil  water  content  near  the  surface  while  leaving  the 
deeper  substrate  relatively  moist.  Jackson  (1973),  CC97, 
and  Carlson  et  al.  (2000)  show  that  this  soil  drying  may 
be  confined  for  a  while  to  a  very  shallow  substrate  layer, 
perhaps  less  than  2  cm  deep.  Decoupling  takes  place 
during  stage-2  (soil  limited)  evaporation  (Idso  et  al. 
1974;  Brutsaert  and  Chen  1995)  in  which  evaporation 
rates  are  limited  by  the  efficiency  of  moisture  diffusion 
from  deeper  to  near-surface  soil  layers.  Under  these  cir- 
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cumstances,  surface  estimates  of  SWC  no  longer  rep¬ 
resent  a  deep  soil  column,  and  current  land  surface  pro¬ 
cess  models  may  have  insufficient  vertical  resolution  to 
capture  this  kind  of  SWC  profile,  so  that  predicted  af¬ 
ternoon  temperatures  may  be  grossly  underestimated. 
Also,  until  recently,  most  atmospheric  prediction  models 
have  simply  employed  climatological  estimates  to  as¬ 
sign  initial  SWC. 

Accurate  initial  vertical  profiles  of  SWC  are  therefore 
necessary  for  making  a  correct  prediction  of  surface 
fluxes.  Hydrological  models  such  as  the  Soil  Hydrology 
Model  (SHM)  used  in  this  study  (Capehart  and  Carlson 
1994;  CC97;  Crawford  et  al.  2000)  can,  in  principle, 
provide  a  better  estimate  of  SWC  profiles  (from  the 
surface  down  to  2  m)  than  can  climatic  data,  using  con¬ 
ventional  meteorological  data.  Moreover,  they  can  cap¬ 
ture  small  variations  in  soil  drying  that  would  otherwise 
produce  serious  errors  in  the  surface  fluxes  if  ignored. 

Fractional  vegetation  cover  describes  the  amount  of 
solar  radiation  incident  on  bare  soil  surfaces  within  a 
specified  area.  For  example,  an  Fr  of  0.75  for  a  particular 
grid  cell  indicates  that  75%  of  that  grid  cell  is  covered 
by  vegetation  and  25%  by  bare  soil,  which  information 
land  surface  schemes  then  apply  to  the  partitioning  of 
surface  processes  such  as  net  radiation  (canopy  vs  soil), 
sensible  heat  flux,  evaporation,  and  transpiration.  Over 
dense  vegetation,  the  latent  heat  flux  is  largely  influ¬ 
enced  by  the  plant  (rather  than  the  underlying  soil  sur¬ 
face)  and  responds  more  sensitively  to  the  SWC  in  the 
root  zone  than  at  the  soil  surface.  Decoupling  events 
are  thus  expected  to  occur  in  low  to  moderately  vege¬ 
tated  areas,  although  they  may  nevertheless  take  place 
at  subgrid  scales  in  bare  patches  between  plants. 

Fractional  vegetation  cover  is  a  parameter  utilized  in 
most  land  surface  schemes;  however,  it  is  almost  always 
specified  using  climatological  data.  The  advent  and  im¬ 
provement  of  remotely  sensed  vegetation  indices  and 
the  close  relationship  between  Normalized  Difference 
Vegetation  Index  (NDVI)  data  and  Fr  have  enabled  more 
accurate  estimations  of  Fr  to  be  produced  at  varying 
temporal  and  spatial  scales  (Choudhury  et  al.  1994;  Gil¬ 
lies  et  al.  1997).  These  estimates  of  Fr  must  also  be 
considered  in  prediction  schemes  in  order  to  capture 
adequately  the  process  of  decoupling. 

This  paper  addresses  two  questions:  1)  under  what 
circumstances  will  decoupling  occur  and  2)  whether 
specifying  more  detailed  and  vertically  resolved  SWC 
and  better  estimates  of  Fr  can  capture  the  decoupling 
process  in  a  coupled  mesoscale  model  [The  Fifth-Gen¬ 
eration  Pennyslvania  State  University  (PSU)“National 
Center  for  Atmospheric  Research  (NCAR)  Mesoscale 
Model  (MM5)]  with  a  land  surface  component  [Simu¬ 
lator  for  Hydrology  and  Energy  Exchange  at  Land  Sur¬ 
face  (SHEELS)].  Accordingly,  we  will  first  discuss  the 
methodology  of  how  initial  SWC  and  Fr  fields  are  gen¬ 
erated.  A  more  detailed  description  of  input  fields  and 
soil  physics  of  MM5/SHEELS  is  included  in  section  2. 
Sections  3  and  4  present  the  output  from  MM5/ 


SHEELS,  with  particular  attention  to  the  impact  of  SWC 
and  Fr  on  the  simulations.  A  validation  of  air  temper¬ 
ature  simulations  is  described  in  section  5.  Section  6 
discusses  these  results  and  complicating  issues  and  as¬ 
sesses  the  value  of  including  these  fields  and  finely  re¬ 
solved  soil  layers  in  land  surface  schemes. 

2,  Generation  of  input  fields 

a.  Soil  water  content 

SHM  was  specifically  designed  to  provide  initial  soil 
water  content  fields  to  hydrological  and  meteorological 
forecast  models  using  only  conventional  land  use  and 
meteorological  data  (Smith  et  al.  1994;  Capehart  1996). 
Its  advantages  over  using  remotely  sensed  or  climato¬ 
logical  data  for  these  mesoscale  models  make  it  desir¬ 
able  for  producing  regional-scale  analyses  of  soil  water 
content  on  a  daily  basis.  More  important,  SHM  can 
provide  a  highly  resolved  vertical  profile  of  soil  mois¬ 
ture  for  any  combination  of  soil  layers.  Although  it  is 
a  one-dimensional  model,  SHM  is  used  here  to  generate 
SWC  fields  at  36-km  grid  spacing  over  most  of  the 
United  States,  from  which  MM5/SHEELS  can  be  ini¬ 
tialized.  Previously,  a  20-layer  version  of  SHM  with  10- 
cm  vertical  spacing  was  used  to  initialize  soil  moisture 
in  MM5/Biosphere-Atmosphere  Transfer  Scheme 
(BATS;  W.  J.  Capehart  1997,  personal  communication; 
Lapenta  et  al.  1998)  and  demonstrated  improvement  in 
predicting  surface  fluxes  over  those  employing  clima¬ 
tological  (default)  SWC  values. 

Still,  this  vertical  spacing  (10  cm)  for  soil  moisture 
may  produce  large  errors  in  surface  fluxes,  which  de¬ 
pend  on  knowledge  of  the  surface  soil  temperature  in 
sparsely  vegetated  areas.  Accordingly,  a  new  higher- 
resolution  9-layer  soil  structure  was  specified  in  SHM 
in  which  the  topmost  layer  is  2  cm  deep  (0-2  cm). 
Subsequent  layers  are  staggered  over  different  substrate 
depths:  3  (2-5),  5  (5-10),  15  (10-25),  25  (25-50),  25 
(50-75),  25  (75-100),  50  (100-150),  and  50  cm  (150- 
200  cm  depth),  for  a  total  of  2  m  in  depth.  This  new 
setup  has  three  layers  in  the  top  1 0  cm  of  soil,  as  opposed 
to  a  single  10-cm  layer.  The  loss  of  resolution  in  the 
middle  and  lower  soil  layers  is  not  thought  to  be  a 
problem  because  of  the  more  uniform  nature  of  SWC 
in  those  deeper  layers.  Furthermore,  tests  on  SHM  have 
shown  that  subgrid-scale  variability  is  smallest  in  deeper 
soil  layers  (Smith  et  al.  1994). 

b.  Fractional  vegetation  cover 

NDVI  is  defined  as 

NDVI  =  (a,R  -  avisV(«iR  +  «vis).  (U 

where  aiR  and  avis  ^^e  the  infrared  and  visible  reflec¬ 
tances,  respectively.  A  square  root  relationship  between 
NDVI  and  Fr  was  first  derived  by  Choudhury  et  al. 
(1994)  and  Gillies  and  Carlson  (1995),  and  subsequently 
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used  by  Gillies  et  al  (1997)  and  Kustas  and  Norman 
(1999).  From  this  result,  it  is  now  possible  to  derive 
reasonable  estimates  of  Fr  directly  from  NDVI.  Carlson 
and  Ripley  (1997)  show  that  calculating  Fr  from  NDVI 
reduces  the  need  for  making  an  atmospheric  correction 
to  the  remote  sensing  data.  For  this  work,  Fr  fields  were 
determined  from  the  Earth  Resources  Observing  Sat¬ 
ellite  (EROS)  Data  Center’s  15-day  Advanced  Very 
High  Resolution  Radiometer  composites  of  NDVI  for 
the  period  1990-98  using  the  method  described  by  Gil¬ 
lies  and  Carlson  (1995).* 

An  example  of  a  final  product  of  NDVI-derived  frac¬ 
tional  vegetation  cover  at  36-km  grid  spacing  over  the 
central  United  States  in  July  can  be  seen  in  Fig.  la.  In 
contrast,  the  climatological  analysis  for  this  period,  as 
generated  by  MM5  and  based  solely  on  subsurface  tem¬ 
perature  and  date,  illustrates  the  lower  horizontal  vari¬ 
ability  of  default  estimates  of  Fr  (Fig.  lb). 

c.  SHEELS  soil  physics 

The  PSU-NCAR  MM5  (Anthes  and  Warner  1978) 
was  chosen  as  the  forecast  model  for  this  study  because 
of  its  wide  use  throughout  the  community  and  its  ease 
of  coupling  with  land  surface  schemes  such  as  BATS 
(Dickinson  et  al.  1986).  Recent  modifications  to  the 
BATS  land  surface  scheme  have  resulted  in  SHEELS 
(Lapenta  et  al.  1998).  Of  crucial  importance  in  this  mod¬ 
el  is  the  inclusion  of  highly  resolved  soil  layers.  Rather 
than  the  0-10-cm  top  layer  of  BATS,  SHEELS  enables 
the  soil  properties,  including  soil  moisture,  to  be  spec¬ 
ified  in  1-cm  layers.  Accordingly,  we  specify  five  soil 
layers  of  equal  depth  in  the  upper  10  cm:  0-2,  2-4,  4-6, 
6-8,  and  8-10  cm.  This  arrangement  can  capture  the 
details  of  rapid  soil  drying  and  the  decoupling  process. 
As  in  BATS,  the  root-zone  depths  are  variable  and  based 
upon  the  land  cover  type. 

The  soil  physics  in  SHEELS  is  nearly  identical  to  that 
in  BATS  (and  SHM),  with  the  exception  of  the  evap¬ 
oration  scheme.  The  temporal  change  in  soil  water  con¬ 
tent  in  each  of  the  soil  sublayers  is  determined  by  par¬ 
titioning  infiltration,  evaporation,  transpiration,  diffu¬ 
sion,  and  gravitational  drainage  (and  surface  runoff  and 
ponded  water  as  residuals).  The  change  in  depth  of  water 
in  soil  layer  with  time  t  can  be  expressed  as 

ddjbt  =  1;  -  E,  -  r,,  (2) 

where  and  T,-  are  the  amounts  of  infiltration,  evap¬ 
oration,  and  transpiration  attributed  to  layer  i  according 
to  weighting  functions.  The  proportion  of  infiltration  for 
each  layer  is  based  on  an  inverted  triangular  weight 
extending  from  the  surface  to  a  depth  determined  by 
land  cover  type.  The  weighting  functions  for  evapora- 


'  A  database  of  all  the  vegetation  cover  fields  has  been  produced 
and  at  the  time  of  writing  was  available  online  at  http://www.essc. 
psu.edu/frac_veg. 


tion  and  transpiration  are  similar  but  are  modified  in 
proportion  to  the  SWC  within  each  layer.  This  approach 
allows  for  more  evapotranspiration  to  be  extracted  from 
wet  layers  than  from  drier  layers,  relative  to  the  weight¬ 
ing  functions,  which  decrease  from  the  surface  down¬ 
ward.  For  evaporation,  the  depth  of  the  triangle  extends 
down  to  the  same  land  cover-dependent  depth  as  for 
infiltration,  while  the  transpiration  depth  extends  down 
through  the  root  zone.  Once  the  water  depths  have  been 
calculated  at  each  time  step,  the  volumetric  water  con¬ 
tents  0,  are  updated. 

At  each  model  time  step  and  at  each  grid  cell,  the 
vertical  fluxes  of  water  between  each  layer  are  calcu¬ 
lated  using  Darcy’s  law: 

=  -kdhfhz,  (3) 

where  =  vertical  water  flux,  k  =  =  hy¬ 

draulic  conductivity,  and  h  =  /z(©)  =  total  hydraulic 
potential  =  if/ +  z  =  ^  [where  if/  =  if/{S)  = 

hydraulic  matric  potential  and  z  =  gravitational  poten¬ 
tial,  b  =  beta  parameter  (reciprocal  of  the  pore  size 
parameter  A),  and  O  =  relative  soil  water  content  = 
(0.  -  0,)/(0,  -  ©r)]'  The  subscripts  5  and  r  refer  to 
saturated  and  residual  (wilting  point)  values,  which  are 
a  function  of  soil  type. 

These  formulations  are  based  on  the  empirical  pa- 
rameterizations  of  soil  properties  by  Clapp  and  Horn- 
berger  (1978),  in  which  the  exponential  value  b  is  rec¬ 
ommended  for  each  designated  soil  type.  Therefore, 

q^  =  -kdif/fdz  -  k 

-  -  /:,(©)^^"^  (4) 

Applying  mass  continuity  yields  Richards’s  equation: 

dd/dt  =  -dqjdz  =  (Bldzykdif/ldz)  +  dk/dz,  (5) 

and  expressing  terms  as  functions  of  0,  we  obtain 

do/dt  =  (d/dz)[D{e)de/dz]  +  (dkide)defdz,  (6) 

where  D{0)  =  kdif/iBO  =  the  diffusion  coefficient.  The 
nonlinear  dependence  of  SWC  fluxes  on  the  hydraulic 
conductivity  and  matric  potential  (or  SWC  itself),  and 
therefore  the  beta  parameter,  is  the  cause  of  rapid  soil 
drying  and  will  be  discussed  further. 

Isothermal  vapor  flux  and  thermally  driven  liquid  and 
vapor  fluxes  occur  near  the  soil  surface.  It  should  be 
noted  that  neither  the  SHM  nor  SHEELS  accounts  for 
these  processes,  and  would  require  the  addition  of  a  heat 
diffusion  equation  to  the  liquid  and  vapor  conservation 
equations,  thus  greatly  increasing  computational  de¬ 
mand  and  complexity  to  a  coupled  model  such  as  MM5/ 
SHEELS.  Studies  have  shown  that  in  the  absence  of 
very  moist  or  desiccated  soil  conditions,  neglect  of  these 
thermally  driven  flows  has  only  a  1%  impact  on  surface 
evaporation,  and  isothermal  liquid  flux  is  the  limiting 
process  on  evaporation  (Milly  1982,  1984;  Saravana- 
pavan  and  Salvucci  2000).  Furthermore,  BATS  and 
SHEELS  resolve  only  a  surface  layer  of  20  cm  and  deep 
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Fig.  L  Fractional  vegetation  cover  over  the  central  United  States  at  36-km  grid  spacing  for 
5-19  Jul  1996,  derived  from  (a)  2-week  composite  NDVI  measurements  and  (b)  climatological 
estimates  generated  by  default  in  MM5.  Values  range  from  0.0  (0%)  to  LO  (100%)  coverage, 
with  brighter  shades  indicating  higher  values. 


soil  layer  temperature  using  a  force-restore  approxi¬ 
mation,  which  is  not  nearly  enough  to  estimate  realistic 
thermal  gradients  and  fluxes.  Results  will  show  that  the 
soil  maintains  a  middle  range  of  SWC  throughout  our 
simulations  and  does  not  reach  saturated  or  desiccated 
conditions.  For  these  reasons,  the  omission  of  thermally 
driven  water  flow  in  these  simulations  is  not  of  primary 
concern  and  is  not  expected  to  produce  first-order  errors. 


A  more  complex  and  robust  version  of  SHEELS  will 
include  a  diffusion-based  soil  scheme  in  which  tem¬ 
peratures  will  be  evaluated  at  each  sublayer,  and  the 
significance  of  these  fluxes  can  be  examined  under  a 
wide  range  of  conditions. 

The  hysteretic  nature  of  wetting  and  drying  fronts 
observed  in  soils  is  not  easily  modeled,  even  by  explicit 
soil  hydrology  models.  Although  hysteresis  does  have 
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Table  L  MM5/SHEELS  simulations  and  the  associated  input  fields 
initialized  in  each  of  the  four  simulations  reported  in  the  paper. 


Simula¬ 

tion 

SHEELS 

physics 

GLCC 

land 

cover 

Fractional 
STATSGO  vegetation 
soil  texture  cover 

SHM 

soil 

moisture 

1 

No 

No 

No 

No 

No 

2 

Yes 

Yes 

Yes 

No 

No 

3 

Yes 

Yes 

Yes 

Yes 

No 

4 

Yes 

Yes 

Yes 

Yes 

Yes 

an  impact  on  soil  water  retention  (j/^0  relationship) 
during  sequential  drying  and  wetting  events,  the  36-h 
simulation  of  a  portion  of  one  drying  event  (or  a  single 
scanning  curve)  performed  here  should  not  be  a  concern. 

d.  MM5/SHEELS  simulations 

MM5/SHEELS  was  run  on  three  nested  domains,  a 
coarse  outer  grid  of  36-km  spacing  (Figs.  la,b),  an  in¬ 
terior  grid  of  12  km,  and  a  finer  mesh  of  4-km  grid 
spacing.  The  4-km  grid  corresponds  to  the  Atmospheric 
Radiation  Measurement  Program  Cloud  and  Radiation 
Test  Bed  (ARM  CART)  Southern  Great  Plains  (SGP) 
site  covering  central  Kansas  and  Oklahoma,  a  region 
characterized  by  a  large  vegetation  gradient  from  east 
to  west.  As  such,  the  domain  is  optimal  for  this  study, 
because  the  large  and  heterogeneous  spatial  variations 
in  Fr  and  SWC  provide  a  suitably  large  range  of  con¬ 
ditions  for  capturing  the  decoupling  effect. 

Dates  were  chosen  from  a  period  of  abundant  sun¬ 
shine  and  light  winds  immediately  following  a  moderate 
precipitation  event.  Initializing  the  model  in  these  con¬ 
ditions  provides  sufficient  time  for  MM5/SHEELS  to 
simulate  decoupling  from  a  moist,  but  not  saturated,  soil 
profile.  Accordingly,  1200  UTC  16  July  1996  was  cho¬ 
sen  as  the  start  time  of  36-h  simulations  spanning  2  days 
of  strong  solar  heating  and  soil  drying.  National  Centers 
for  Environmental  Prediction  reanalysis  data  were  used 
as  the  atmospheric  boundary  conditions  and  are  not  a 
major  concern  in  our  analyses. 

Four  MM5/SHEELS  simulations  for  this  36-h  period 
were  made:  1)  MM5  coupled  with  BATS  using  all  cli¬ 
matological  (default)  parameters,  2)  MM5  coupled  with 
SHEELS  using  default  vegetation  cover  and  soil  mois¬ 
ture  profiles,  3)  MM5  coupled  with  SHEELS  using  the 
derived  Fr  and  default  soil  moisture  profiles,  and  4) 
MM5  coupled  with  SHEELS  using  the  derived  Fr  and 
derived  soil  moisture  profiles  from  SHM.  These  simu¬ 
lations  and  input  fields  are  summarized  in  Table  1. 

MM5/SHEELS  divides  its  flux  computation  into  a 
fraction  of  each  grid  cell  covered  by  vegetation  and  a 
fraction  covered  by  bare  soil,  then  averages  the  two  to 
produce  a  representative  value  for  each  grid  cell  (Smith 
et  al.  1994).  Fractional  vegetation  cover  was  initialized 
for  these  simulations  using  the  appropriate  2-week 
NDVI  composite  and  processed  as  described  in  section 
2b.  The  derived  Fr  data  were  aggregated  to  the  36-km 


MM5  grid  using  a  spatial  average  of  1296  1-km  cells 
for  each  MM5  cell. 

SHM  and  MM5  use  the  State  Soil  Geographic  Da¬ 
tabase  (STATSGO)  dataset  of  soil  types  (Miller  et  al. 
1994),  where  1  of  12  soil  types  is  assigned  to  each  1 
km  X  1  km  cell  covering  the  entire  United  States.  Each 
of  the  12  soil  types  is  assigned  values  for  hydraulic 
properties  such  as  minimum  soil  suction,  saturated  con¬ 
ductivity,  porosity,  wilting  point,  and  the  beta  parameter, 
which  the  models  then  use  to  calculate  the  vertical  fluxes 
of  soil  moisture  (section  2c).  Land  cover  data  were  gen¬ 
erated  from  the  Global  Land  Cover  Characteristic 
(GLCC)  project  operated  by  the  EROS  Data  Center  in 
1997  (http://edcdaac.usgs.gov/glcc/glcc.html).  Each 
grid  cell  is  assigned  1  of  18  land  cover  classes  that 
correspond  to  biophysical  parameter  estimates,  accord¬ 
ing  to  a  lookup  table,  such  as  leaf  area  density,  rough¬ 
ness  length,  and  albedo  (BATS;  Dickinson  et  al.  1986). 

Initial  fields  of  soil  type  (STATSGO)  and  land  cover 
(GLCC)  were  used  in  both  the  SHM  and  MM5  simu¬ 
lations  at  36-km  grid  spacing.  Each  was  aggregated  from 
1  to  36  km  by  assigning  the  most  frequent  class  or  type 
occurring  in  each  36-km  grid  cell.  These  fields,  along 
with  the  36-km  derived  Fr  and  SWC  fields  (simulations 
3  and  4),  were  uniformly  interpolated  by  MM5  down 
to  the  grid  spacing  of  the  12-  and  4-km  nests  upon 
initialization.  Although  the  land  cover,  soil  type,  and  Fr 
data  are  available  at  1-km  grid  spacing,  we  felt  it  crucial 
to  keep  consistency  between  SHM  and  MM5.  This  en¬ 
sured  that  SWC  estimated  by  SHM  for  a  particular  grid 
cell  of,  say,  short  grass,  loam,  and  26%  Fr  would  be 
simulated  in  MM5  under  the  same  conditions. 

In  SHM,  the  Clapp  and  Hornberger  (CH;  1978) 
scheme  is  used,  a  portion  of  which  is  shown  in  the  top 
part  of  Table  2,  wherein  the  12  STATSGO  soil  types 
and  their  associated  soil  parameter  values  are  given. 
BATS  and  SHEELS,  however,  use  the  BATS  table  of 
hydraulic  constants  (bottom  of  Table  2),  which  differs 
from  that  of  CH.  BATS  allows  SWC  to  vary  from  the 
wilting  point  (residual  SWC)  of  the  soil,  which  ranges 
from  0.029  to  0.3575,  to  saturation,  and  CH  allows  SWC 
to  vary  from  0  to  saturation.  In  addition,  the  values  of 
saturation  in  BATS  and  CH  are  not  equal.  Therefore, 
the  range  of  allowable  SWC  may  differ  significantly 
between  models,  and  a  rescaling  of  the  soil  moisture 
output  from  SHM  was  performed  before  it  could  be 
included  in  MM5  simulations: 


(5)  = 

'JmM5  1  ^  1 

'-’s.CH  / 


0. 


v@ 


—  (h) 

.s,BATS  '-^r.BATS 


)  + 


BATS> 


(7) 


where  0,,shm  is  the  soil  water  content  from  SHM,  and 
the  saturation  (subscript  s)  and  residual  values  (subscript 
r)  are  obtained  from  the  BATS  and  CH  soil  parameter- 
izations  (Table  2).  The  SWC  in  the  upper  10-cm  layer 
of  soil,  if  estimated  as  0.235  from  SHM,  would  be  as¬ 
signed  a  value  of  0.261  for  a  silt  loam  soil  upon  ini¬ 
tialization  in  MM5/SHEELS  using  this  technique.  The 
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Table  2.  Upper  (saliirated)  and  lower  (residual)  limiting  values  of 
soil  water  content  (O)  and  the  beta  parameters  for  the  12  STATSGO 
soil  types,  as  specified  by  Clapp  and  Hornberger  (1978;  CH)  and 
BATS  (Dickinson  et  al.  1986). 


Soil  type 

SWC,,  (0,3 

SWC„„„  (0,) 

h 

CH 

Sand 

0.395 

0 

3.5 

Loamy  sand 

0.41 

0 

4.0 

Sandy  loam 

0.435 

0 

4.5 

Silt  loam 

0.485 

0 

5.0 

Sill 

0.485 

0 

5.5 

Loam 

0.451 

0 

6.0 

Sandy  clay  loam 

0.42 

0 

6.8 

Silly  clay  loam 

0.477 

0 

7.6 

Clay  loam 

0.476 

0 

8.4 

Sandy  clay 

0.426 

0 

9.2 

Silly  clay 

0.492 

0 

10.0 

Clay 

0.482 

0 

10.8 

BATS 

Sand 

0.33 

0.029 

3.5 

Loamy  sand 

0.36 

0.0428 

4.0 

Sandy  loam 

0.39 

0.0588 

4.5 

Silt  loam 

0.42 

0.1117 

5.0 

Silt 

0.45 

0.1350 

5.5 

Loam 

0.48 

0.1592 

6.0 

Sandy  clay  loam 

0.51 

0.1926 

6.8 

Silty  clay  loam 

0.54 

0.2260 

7.6 

Clay  loam 

0.57 

0.2592 

8.4 

Sandy  clay 

0.60 

0.2922 

9.2 

Silty  clay 

0.63 

0.3250 

10.0 

Clay 

0.66 

0.3575 

10.8 

ramifications  of  this  rescaling  between  models  will  be 
discussed  in  section  6. 

3.  Simulation  results:  SWC 

a.  Soil  moisture  initialization 

Although  MM5/SHEELS  resolves  soil  properties  in 
2-cm  layers,  SWC  must  be  initialized  in  the  model  in 
the  original  BATS  layers  of  an  upper  (0-10  cm),  root 
zone,  and  deep  soil  layer  down  through  2  m.  The  9-layer 
SWC  data  from  SHM  were  aggregated  into  these  3  lay¬ 
ers  for  use  in  simulation  4,  and  the  upper  10-cm  layer 
at  initial  time  is  shown  in  Fig.  2a  for  the  4-km  nested 
domain.  This  aggregation  of  SHM  data  and  loss  of  ver¬ 
tical  resolution  upon  initialization  will  be  discussed  in 
section  6. 

In  comparison,  the  default  0-10-cm  SWC  initialized 
by  MM5  in  simulations  1,  2,  and  3  is  shown  in  Fig.  2b. 
The  default  SWC  is  purely  a  function  of  land  cover  type 
and  date  and  is  initialized  from  a  lookup  table  in  MM5/ 
SHEELS.  This  default  estimate  gives  a  single  value  of 
SWC  that  is  valid  for  the  entire  soil  column  at  each  grid 
cell,  from  the  surface  down  through  2  m,  and  therefore 
lacks  any  vertical  variability,  the  implications  of  which 
will  be  discussed  in  section  6.  It  is  apparent  that  the 
spatial  variability  of  soil  moisture  from  SHM  (Fig.  2a) 
is  much  higher  than  the  default  field  generated  for  this 
date.  The  sharp  gradients  of  SWC  in  both  figures  rep¬ 


resent  the  soil  type  boundaries,  corresponding  to  the 
gradients  in  hydraulic  properties  governed  by  soil  type. 

b.  Near-surface  soil  drying 

Soil  moisture  simulations  show  that  some  locations 
dry  out  quickly  in  the  upper  part  of  the  five  sublayers 
while  others  remain  wet.  Indeed,  a  large  variability  of 
drying  rates  can  be  seen  throughout  the  ARM  CART 
region  in  these  simulations.  Figure  3  shows  two  simu¬ 
lations  of  the  upper  five  sublayers  of  SWC  using  derived 
SWC,  labeled  SHM,  in  one  case  and  the  default  SWC, 
labeled  FR,  in  another  at  Lamont,  Oklahoma.  The  rep¬ 
resentative  soil  type  here  is  listed  as  silt  loam,  and  the 
vegetation  fraction  is  0.187.  The  lowest  (driest)  curves 
in  each  figure  represent  the  uppermost  (0-2  cm)  soil 
layer.  Through  the  first  24  h  of  the  simulation,  the  de¬ 
fault  SWC  shows  much  greater  signs  of  decoupling,  as 
the  upper  2  cm  dries  out  nearly  0.04  more  than  does 
the  8-10-cm  layer.  On  the  other  hand,  the  derived  SWC 
shows  little  variation  in  drying  between  layers  through 
24  h.  Only  late  on  the  second  afternoon  does  the  sim¬ 
ulation  with  the  derived  SWC  begin  to  show  more  sig¬ 
nificant  drying  near  the  surface. 

Figure  4  shows  the  upper  five  soil  layers  at  Coldwater, 
Kansas,  for  the  same  two  simulations.  Coldwater  is  also 
located  in  a  relatively  bare  soil  region,  with  a  vegetation 
cover  of  0.257  and  a  soil  type  of  loam.  What  is  partic¬ 
ularly  noticeable  at  this  location  is  the  great  amount  of 
decoupling  seen  in  using  the  derived  values  of  SWC 
through  the  36-h  period.  The  derived  SWC  shows  a 
decrease  in  the  0-2-cm  layer  of  nearly  0.07  below  that 
of  the  8-10-cm  layer  and  0.03  below  that  of  the  2-4- 
cm  layer. 

El  Reno,  Oklahoma,  is  another  sparsely  vegetated  lo¬ 
cation  with  an  Fr  of  0.247  and  silt  loam  soil.  The  soil 
moisture  curves  for  the  two  simulations  for  El  Reno  are 
shown  in  Fig.  5.  Here,  the  simulation  with  the  derived 
SWC  starts  out  much  wetter  again  but  shows  no  indi¬ 
cation  of  decoupling  at  all  during  the  36-h  simulation. 
The  default  SWC  shows  moderate  decoupling,  similar 
to  that  seen  in  Lamont,  which  has  identical  soil  type 
and  similar  vegetation  cover  properties.  The  noticeable 
soil  drying  differences  of  the  derived  SWC  simulations 
between  Lamont  and  El  Reno,  particularly  on  day  2,  are 
due  to  the  initial  SWC  values.  The  derived  SWC  values 
at  Lamont  in  the  0-10-cm  layer  are  about  0.26;  in  El 
Reno  the  initial  SWC  is  near  0.29. 

c.  Drying  thresholds  and  soil  type 

At  Lamont,  the  derived  SWC  values  lead  to  decou¬ 
pling  on  the  second  afternoon  because,  as  suggested  by 
CC97,  the  SWC  reaches  a  threshold  level  for  rapid  dry¬ 
ing  in  this  soil  type,  whereas  in  El  Reno  the  derived 
SWC  remains  too  high  for  rapid  drying  to  occur.  Upon 
investigation  we  find  that  the  degree  of  surface  drying 
depends  crucially  on  the  initial  SWC  values.  Figures 
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Fig.  2.  Soil  water  content  (0)  of  the  O-lO-cm  layer  initialized  in  MM5/SHEELS  at  1200  UTC 
16  Jul  1996,  derived  from  (a)  SHM  soil  moisture  profiles  and  (b)  climatological  estimates  gen¬ 
erated  by  default  in  MM5.  This  4-km  nested  grid  corresponds  to  that  of  the  ARM  CART  project 
over  central  Oklahoma  and  Kansas. 
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Fig.  3.  Soil  water  content  (0)  of  the  upper  five  sublayers  over  the  36-h  period  for  Lament, 
OK.  The  open  circles  (SHM)  represent  the  SHM-derived  SWC  simulation,  the  default  SWC 
simulation  is  represented  by  closed  circles  (FR).  The  0-2-cm  layer  (lowest  curve)  is  the  driest  in 
each,  and  the  S-lO-cm  layer  (highest  curve)  is  the  wet  test. 


3-5  show  that  MM5/SHEELS  simulates  rapid  drying  of 
near-surface  soil  layers  only  when  the  upper  10  cm  of 
soil  moisture  resides  in  a  certain  threshold  range  that  is 
highly  dependent  on  the  soil  type. 

Thresholds  for  different  soil  types  can  be  determined 
by  comparing  drying  rates  and  SWC  values  with  soil 
type  throughout  the  inner  nest.  When  the  rate  of  change 


of  SWC  for  the  0-2-cm  layer  is  significantly  greater 
than  that  of  the  8-10-cm  layer,  the  value  of  SWC  at 
which  this  divergence  in  drying  rate  first  occurs  is  noted; 
this  establishes  an  upper  threshold  range  for  a  given  soil 
type.  A  similar  procedure  determines  the  lower  thresh¬ 
old  limit.  Table  3  gives  the  thresholds  estimated  from 
these  simulations  for  the  six  soil  types  present  in  our 


Fig.  4.  Same  as  for  Fig.  3  but  for  Coldwater,  KS. 
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domain.  Classes  such  as  silty  clay  loam  and  clay  loam 
(soil  types  8  and  9)  tend  to  have  higher  threshold  levels 
due  to  their  higher  saturation  values.  Consequently,  sim¬ 
ulations  using  the  higher  initial  values  of  SWC  (the 
simulations  using  the  derived  SWC)  experience  a  more 
pronounced  decoupling  for  soil  types  8  and  9.  The  sim¬ 
ulation  with  default  soil  moisture  and  derived  Fr  shows 
more  rapid  drying  for  soil  types  1—4  because  the  typical 
default  values  of  SWC  lie  in  these  lower  threshold  rang¬ 
es. 

Table  3  also  lists  the  percent  of  saturation  and  field 
capacity  for  each  soil  type  and  corresponding  threshold 
range.  The  upper  limits  of  these  thresholds  range  from 
39.2%  to  50.2%  of  saturation,  and  all  but  silt  loam 
(49.1%)  have  ranges  that  include  values  of  50%  of  field 
capacity.  This  result  is  in  agreement  with  an  earlier  study 
by  CC97,  who  found  similar  thresholds  for  decoupling 
events  using  a  different  land  surface  scheme  and  hy¬ 
draulic  equations.  Overall,  rapid  drying  is  likely  if  the 
initial  SWC  is  between  13%  and  50%  of  saturation, 
depending  on  the  soil  type  and  soil  scheme  used,  and 


an  SWC  near  50%  of  field  capacity  almost  always  leads 
to  rapid  drying. 

CC97  explain  how  the  soil  physics  scheme  [Eqs. 
(4-7)1  governing  the  diffusion  and  evaporation  of  SWC 
can  simulate  decoupling.  According  to  these  equations, 
hydraulic  conductivity  varies  with  water  content  to  the 
power  of  3  +  lb,  and  matric  potential  varies  with  water 
content  to  the  —b  power,  where  b  is  empirically  derived 
and  assigned  by  soil  type  according  to  CH  (1978).  When 
a  moist  column  of  soil  begins  to  evaporate,  the  hydraulic 
conductivity  k  decreases  slowly  while  the  matric  po¬ 
tential  gradient  dtp/dz  increases  slowly.  For  very  moist 
and  dry  values  of  ft  these  changes  in  k  and  ip  tend  to 
offset  each  other,  keeping  a  constant  upward  flux  of 
moisture  and  preventing  rapid  drying.  However,  as 
evaporation  increases  throughout  the  day,  6  enters  a 
middle  range  of  SWC  where  the  conductivity  tends  to 
decrease  rapidly  (according  to  the  3  +  2b  power),  and 
the  matric  potential  gradient  does  not  increase  rapidly 
enough  to  compensate  (according  to  the  -b  power 
only).  Therefore,  the  upward  flux  of  moisture  from  the 


Table  3.  Threshold  values  (0^,^)  and  corresponding  relative  saturation  and  field  capacity  percentages  of  the  upper  10  cm  of  SWC  for 
each  soil  type  that  produce  decoupled  drying  profiles  of  the  upper  soil  layers. _ 


Soil  texture  class 

0, 

(BATS— Table  2) 

Initial  SWC 
of  drying  events 
(©.,) 

Relative  SWC  saturation  % 
(0,,,  -  0,)/(©.,  -  0.) 

%  of  field  capacity 
(0.,„  -  0.)/O.75(©,  -  ©,) 

1  (Sand) 

3  (Sandy  loam) 

4  (Silt  loam) 

6  (Loam) 

8  (Silty  clay  loam) 

9  (Clay  loam) 

0.33 

0.39 

0.42 

0.48 

0.54 

0.57 

0.10-0.17 

0.170-0.225 

0.175-0.225 

0.235-0.285 

0-27-0.35 

0.30-0.40 

23.6- 46.8 

33.6- 50.2 

20.5- 36.8 

23.6- 39.2 

14.1- 39.5 

13.1- 45.3 

31.5- 62.4 

44.8- 66.9 

27.3-49.1 

31.5- 52.3 

18.8- 52.7 

17.5- 60.4 
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lower  soil  layers  is  retarded  according  to  Eq.  (3).  This 
process  continues  as  the  upper  soil  layer  continues  to 
evaporate  at  a  high  rate  but  is  not  replenished  from 
below,  and  it  therefore  dries  out  considerably  more  than 
the  lower  layers.  The  lower  threshold  limit  occurs  when 
evaporative  demand  decreases  because  of  the  dryness 
of  the  soil  (as  determined  by  the  evaporative  weighting 
scheme). 

Thus,  decoupling  is  not  seen  in  soils  that  are  already 
dry  or  near  saturation.  That  the  exact  threshold  ranges  are 
strictly  a  function  of  the  soil  type  {b)  suggests  that  rapid 
soil  drying  is  sensitive  to  the  choice  of  parameterization. 
Because  the  surface  drying  profoundly  affects  the  surface 
temperature  and  moisture,  differing  threshold  ranges  may 
produce  significant  differences  in  sensible  and  latent  heat 
fluxes.  CC97  show  that  a  change  in  the  h  parameter  in 
the  soil  water  flux  equation  (CH)  of  only  one  standard 
deviation  within  a  single  soil  type  results  in  differences  in 
the  drying  rates  of  the  soil  layers  that  are  almost  as  sig¬ 
nificant  as  those  between  soil  types.  Therefore,  it  is  not 
so  much  the  actual  threshold  values  that  are  of  importance 
in  these  results.  What  these  results  do  show  is  the  close 
dependence  of  initial  SWC  and  drying  rates  simulated  in 
land  surface  schemes  on  the  parameterization  of  soil  type 
and  on  the  choice  of  soil  parameters. 

4.  Simulation  results:  Fr 

a.  Overall  fractional  vegetation  cover  fields 

The  NDVI-derived  Fr  field  over  the  ARM  CART  (4 
km)  domain  for  the  period  5-19  July  1996,  shown  in 
Fig.  6a,  exhibits  great  spatial  detail  in  capturing  the  sharp 
east-west  gradient  of  Fr.  Figure  6b  shows  the  default  Fr 
over  the  same  region  on  16  July  1996.  Default  Fr  is 
initialized  in  MM5/SHEELS  from  a  lookup  table  and  is 
a  function  of  subsurface  soil  temperature,  land  cover  type, 
and  a  seasonal  range  of  vegetation  cover,  as  in  the  original 
BATS  (Dickinson  et  al.  1986).  The  homogeneous  nature 
of  Fr,  which  is  mostly  above  0.8  because  of  the  uniformly 
high  temperatures  of  July,  fails  to  capture  the  actual  spa¬ 
tial  variation  in  vegetation  that  exists  over  the  western 
half  of  the  central  Great  Plains  region. 

h.  Impacts  of  Fr  on  MM5/SHEELS  simulations: 

Temporal  evolution  of  soil  moisture 

Figures  7a,b  show  the  0-2-cm  soil  water  content  for 
simulations  3  and  2,  respectively,  along  the  36.17°N 
latitude  stripe  over  the  36-h  period.  The  simulations  are 
identical  except  for  the  Fr  field  initialized  in  each  (de¬ 
rived,  Fig.  7a,  or  default.  Fig.  7b).  In  both  figures,  sharp 
horizontal  gradients  in  SWC  result  from  the  imposition 
of  soil  type  gradients  across  the  latitude  stripe.  With 
default  values  of  Fr,  the  SWC  remains  high,  and  rapid 
drying  is  absent.  With  derived  Fr,  a  wide  area  of  rapid 
drying  is  evident  in  the  middle  of  the  segment,  and  a 
check  of  the  SWC  profiles  in  the  region  confirms  that 


decoupling  does  occur  here.  This  central  area  of  drying 
corresponds  to  a  soil  type  of  silt  loam  and  to  an  initial 
SWC  (identical  in  each  run)  within  the  corresponding 
threshold  range  for  decoupling.  The  derived  Fr  gives 
less  than  0.3  (30%)  vegetation  cover  (Fig.  8)  for  this 
central  region  and  so  results  in  significant  evaporation 
from  a  largely  bare  soil  surface.  The  default  vegetation 
cover  for  this  area  remains  above  0.8  (80%);  the  upper 
soil  layers,  therefore,  do  not  dry  out  as  quickly,  despite 
the  initial  SWC  being  within  the  drying  threshold  range. 

On  the  eastern  edge  of  the  domain,  the  soil  type  is 
also  silt  loam,  and  the  initial  SWC  is  within  the  threshold 
range  for  rapid  drying,  as  in  the  central  region.  However, 
drying  is  much  less  evident  and  much  slower  in  this 
area  with  derived  Fr  than  in  the  central  portion,  with 
only  a  slight  decoupling  (not  shown)  of  the  soil.  This 
result  is  due  to  the  higher  values  of  Fr  (above  0.60) 
along  the  eastern  edge  of  the  cross  section,  which  limits 
the  amount  of  bare-soil  evaporation.  The  remainder  of 
the  cross  section  experiences  little  drying,  because  the 
initial  SWC  values  of  the  simulations  are  not  within  the 
threshold  ranges  for  decoupling. 

c.  Consequences  of  soil  drying  due  to  vegetation 

cover  variations 

Soil  drying  manife.sts  itself  in  the  surface  fluxes  and 
the  surface  temperatures.  Figure  9  shows  the  sensible 
heat  fluxes  across  the  36. 1 T  latitude  stripe  at  1 900  UTC 
on  17  July  1996,  approximately  midafternoon  on  the 
second  day  of  the  simulation  when  considerable  drying 
has  already  occurred.  The  simulation  with  derived  Fr 
shows  generally  larger  sensible  heat  fluxes  than  that 
with  the  default  Fr,  the  greatest  differences  (greater  than 
a  factor  of  2)  occurring  in  the  central  region  where 
decoupling  is  most  pronounced.  Both  simulations  show 
nearly  equal  fluxes  along  the  eastern  border,  where  de¬ 
rived  and  default  Fr  are  both  large.  The  latent  heat  fluxes 
(not  shown)  exhibit  the  opposite  behavior,  as  expected. 

Surface  skin  temperatures,  simulated  for  the  4-km 
domain  at  late  afternoon  on  the  second  day  of  the  sim¬ 
ulations,  are  shown  in  Figs.  10a  (derived  Fr)  and  10b 
(default  Fr).  Derived  Fr  (Fig.  10a)  produces  markedly 
higher  surface  temperatures  than  the  simulation  with 
climatological  Fr  (Fig.  10b),  with  the  exception  of  the 
eastern  and  southeastern  edges.  The  temperature  pattern 
inversely  corresponds  with  the  distribution  of  Fr  (Figs. 
6a, b),  and  therefore  reflects  the  percentage  of  bare  soil 
present  in  each  grid  cell.  The  variability  within  the  west¬ 
ern  half  of  the  domain  in  each  simulation  reflects  the 
happenstance  distribution  of  the  drying  threshold,  which 
is  a  function  of  the  soil  type. 

5.  Validation  of  MM5/SHEELS  simulations 

Air  temperature  measurements  from  the  ARM  CART 
domain  were  examined  on  a  point-by-point  basis  and 
compared  with  simulated  values  for  those  locations.  Fig- 
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Fig.  7.  SWC  along  the  36.17°N  latitude  stripe  of  the  4-km  domain,  showing  0-2-cm  SWC  as 
a  function  of  time  (increasing  upward)  and  distance  for  the  (a)  NDVI-derived  Fr  and  (b)  default 
Fr  simulations  during  the  36-h  period.  Note  the  rapid  drying  evident  in  the  central  portion  of  the 
cross  section  for  the  derived  Fr  simulation. 


ure  1 1  shows  the  temperature  curves  of  the  four  sim¬ 
ulations  and  the  corresponding  observations  over  the 
36-h  period  for  Lamont.  Simulation  3  clearly  yields  the 
closest  agreement,  particularly  in  the  maximum  values. 
The  calculated  rmses  of  the  four  simulations  are  as  fol¬ 
lows:  1  =  2.562,  2  =  2.320,  3  =  1.729,  and  4  =  3.245. 


All  four  simulations  yield  minimum  temperatures  above 
the  observed,  with  those  using  derived  SWC  being  clos¬ 
est,  but  not  within  2'^C  of  the  observed  minimums. 

Soil  moisture  profiles  from  simulations  with  default 
versus  derived  SWC  are  presented  for  this  station  in 
Fig.  3.  The  former  shows  decoupling  at  this  location 
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Fig.  8.  Fractional  vegetation  cover  along  the  36.17°N  latitude  stripe.  The  open  circles  (FR) 
represent  derived  (NDVI  based)  vegetation  cover;  the  closed  circles  (SH)  represent  default  veg¬ 
etation  cover.  Here,  Fr  is  constant  over  the  36-h  period  because  the  2-week  composites  of  NDVI 
data  and  slow  climatological  variability  of  default  estimates  yield  consistent  Fr  throughout. 


throughout  the  36-h  simulation,  and  the  latter  only  be¬ 
gins  to  dry  out  at  the  very  end  of  the  period.  This  de¬ 
coupling  of  the  default  SWC  results  in  a  large  increase 
in  sensible  heat  flux  on  both  afternoons  and,  therefore, 
higher  maximum  temperatures.  The  derived  SWC,  much 
wetter  down  through  the  root  zone,  yields  large  evap¬ 
oration  and  transpiration  rates  on  day  1 .  When  the  upper 


soil  starts  to  dry  out  late  on  the  second  day,  however, 
the  latent  heat  flux  begins  to  diminish  with  time,  and 
the  sensible  heat  flux  increases  (not  shown).  It  is  likely 
that  this  drying  would  have  led  to  decoupling  on  the 
third  day,  thereby  reducing  the  forecast  errors  for  the 
simulation  with  derived  SWC. 

Measurements  of  SWC  at  2.5  cm  (0“5-cm  SWC  av- 


Fig.  9.  Sensible  heat  flux  (W  m  across  the  36.17°N  latitude  stripe  at  1900  UTC  18  Jul  1996. 
Open  circles  (FR)  represent  the  derived  (NDVI-based)  vegetation  cover  simulation;  the  closed 
circles  (SH)  represent  the  default  vegetation  cover  simulation. 
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Fig.  11.  Simulated  (simulations  1-4)  and  measured  (OBS)  air  temperature  for  Lament,  OK, 
over  the  36-h  period.  Refer  to  Table  1  for  the  input  parameters  used  in  simulations  1  (here  labeled 
“MM5”),  2  (“SHE”),  3  (“FR”),  and  4  (“SHM”). 


erages  measured  in  plastic  bags)  from  this  site  show  a 
slow  decrease  from  0.29-  to  0.250,  followed  by  a  more 
rapid  drying  trend  beginning  in  the  0.20-0.23  range  by 
the  second  afternoon.  This  observed  drying  may  reflect 
a  threshold  for  this  soil  type,  consistent  with  the  estimate 
of  0.175-0.225  derived  from  the  simulations.  According 
to  these  measurements,  the  simulation  with  default  SWC 
provides  a  better  temperature  simulation  because  its  val¬ 
ues  happen  to  coincide  with  the  decoupling  threshold. 

The  12  other  ARM  CART  facilities  display  almost 
identical  patterns  to  that  of  Lamont.  The  simulation  with 
default  SWC  and  derived  Fr  yields  good  results  for  the 
afternoon  temperatures;  those  with  the  derived  SWC 
from  the  SHM  were  the  least  accurate,  with  errors  as 
large  as  6°C.  The  fact  that  these  two  runs  differ  only  in 
the  initial  soil  moisture  and  yield  such  dissimilar  results 
underscores  the  importance  of  the  SWC  parameter  in 
land  surface  models.  Even  when  simulations  with  the 
derived  SWC  exhibit  decoupling,  which  occurs  only  on 
day  2  in  most  cases,  temperatures  are  unable  to  reach 
those  higher  values  reported  in  the  observations  because 
of  the  high  evaporation  and  low  temperatures  on  day  1. 
Simulations  1  and  2,  although  they  may  contain  the 
more  accurate  (default)  surface  soil  moisture,  are  greatly 
affected  by  the  high  values  of  default  Fr  (>0.8)  and  are 
consistently  r-2°C  cooler  than  those  with  the  derived 
Fr. 

Overall,  there  was  little  consistent  evidence  of  agree¬ 
ment  between  any  one  simulation  and  the  measured  soil 
moisture,  dewpoint  temperatures,  or  fluxes.  This  result 
may  be  due  to  the  point  versus  area  mismatch  in  scales, 
particularly  with  surface  fluxes,  which  tend  to  have  great 


fluctuations  at  a  point  but  smooth  out  over  larger  areas. 
In  view  of  the  possible  errors  in  the  SWC,  dewpoint, 
and  temperature  measurements  at  these  points,  these  val¬ 
idations  are,  at  best,  tentative,  and  the  only  consistency 
was  seen  in  the  air  temperature  comparisons. 

6.  Discussion  of  results  and  complicating  issues 

a.  Generation  of  soil  moisture  fields 

The  SHM  was  responding  to  a  recent  precipitation 
event,  which  generated  SWC  values  much  higher  than 
the  default  SWC  values  from  MM5/SHEELS.  Because 
BATS  considers  a  higher  range  of  SWC  than  the  CH 
values  (Table  2),  the  soil  moisture  from  SHM  was  scaled 
up  even  more  [Eq.  (1)],  thus  amplifying  the  differences 
between  the  default  and  derived  fields.  Thus,  the  in¬ 
accuracy  of  the  SHM  values,  in  terms  of  resulting  tem¬ 
perature  validations,  may  be  an  artifact  of  the  data  re¬ 
duction  process. 

In  addition  to  the  ambiguity  created  by  scaling,  a 
further  source  of  uncertainty  in  initializing  the  9-layer 
SHM  output  was  that  the  latter  had  to  be  aggregated 
into  three  discrete  layers  for  input  into  MM5/SHEELS, 
the  topmost  being  0-10  cm.  For  this  particular  study, 
not  much  was  lost  in  this  aggregation,  because  the  ini¬ 
tialization  time  is  postprecipitation  when  fairly  uniform 
soil  profiles  were  present  (0.234,  0.235,  and  0.236  at 
Lamont  for  the  top  three  SHM  layers  composing  the 
0-10-cm  layer).  Had  we  been  investigating  an  already 
dry  soil  surface,  such  as  the  one  simulated  by  SHM  for 
18  July,  this  issue  would  have  created  a  much  larger 
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problem,  and  SHEELS  is  being  reformulated  to  account 
for  flexible  soil  layering  at  initialization  as  a  result. 

There  is  also  the  problem  created  by  the  default  SWC 
being  vertically  constant  from  the  surface  through  2  m 
as  generated  by  MM5  for  the  default  SWC  simulations. 
A  look  at  the  FR  simulations  in  Figs.  3-5  shows  a  def¬ 
inite  lack  of  the  diurnal  variability  seen  in  the  SHM 
simulation.  Because  of  the  constant  initial  SWC  profile, 
there  is  not  enough  of  a  gradient  of  SWC  from  the  upper 
layers  to  the  root  zone  to  generate  a  recharge  overnight 
(towards  the  surface)  of  moisture.  The  initial  profile  of 
derived  SWC  is  considerably  more  moist  in  the  root 
zone  and  acts  to  drive  this  recharge  according  to  the 
flow  equations  presented  earlier.  Although  the  default 
profile  simulates  better  temperatures  in  the  end,  this 
result  is  only  due  to  the  initially  drier  near-surface  soil 
layer;  longer  timescales  might  have  yielded  different 
results. 

Regardless  of  which  method  of  determining  SWC  is 
more  accurate,  modeled  or  default,  our  results  demon¬ 
strate  the  importance  of  initial  conditions  on  the  future 
drying  of  the  soil.  Future  studies  linking  hydrologically 
derived  SWC  with  mesoscale  models  should  consider 
using  identical  hydraulic  schemes  to  eliminate  the  need 
for  and  possible  problems  of  rescaling  SWC  data  and 
should  consider  the  use  of  consistent  soil  layers  between 
hydrologic  and  mesoscale  soil  modeling. 

SHM-generated  profiles  do  have  three  distinct  ad¬ 
vantages  over  derived  SWC  profiles  with  regard  to  the 
issues  above:  1)  Modeled  SWC  is  not  climatological 
and  is  driven  by  actual  meteorological  forcings,  2)  There 
is  vertical  variation  in  SWC  profiles  from  SHM,  which 
can  contribute  immediately  to  vertical  moisture  fluxes 
in  a  coupled  model,  and  c)  SHM  can  give  an  indication 
of  when  and  where  decoupling  may  be  likely  by  looking 
at  the  SWC  with  soil  types. 

b.  Generation  of  fractional  vegetation  cover  data 

Clearly,  Fr  estimates  generated  here  are  an  improve¬ 
ment  in  both  time  and  space  over  default  estimates.  As 
evidenced  in  the  results,  an  improvement  in  simulating 
afternoon  temperatures  occurs  when  a  more  detailed 
vegetation  field  is  initialized,  although  the  initial  con¬ 
dition  of  near-surface  soil  moisture  is  often  more  im¬ 
portant  than  the  vegetation  amount  (unless  there  is  near¬ 
ly  complete  vegetation  cover,  in  which  case  transpiration 
and  root-zone  SWC  will  dominate). 

c.  Simulation  of  decoupling 

Results  indicate  that  improved  soil  moisture  data  do 
not  necessarily  yield  more  accurate  simulations  of  land 
surface  processes,  because  of  uncertainties  in  soil  type 
and  vegetation  parameterizations.  Higher  SWC,  such  as 
that  generated  by  the  SHM,  may  lie  above  the  decou¬ 
pling  threshold  range,  which  is  governed  by  the  hy¬ 
draulic  constants  in  SHEELS.  In  that  instance,  the  soil 


may  dry  out  too  slowly  for  the  surface  temperatures  to 
respond  correctly. 

Decoupling  has  not  been  widely  noted  in  mesoscale 
simulations,  suggesting  not  only  that  it  occurs  for  a 
limited  range  of  conditions,  but  that  existing  models 
have  had  insufficient  substrate  resolution  to  describe  it. 
Because  decoupling  is  so  critically  dependent  on  res¬ 
olution,  initial  SWC,  soil  type,  Ff,  and  insolation,  it 
remains  a  challenge  to  predict  where  and  when  it  will 
occur.  This  challenge  is  illustrated  by  Fig.  5.  When  the 
SWC  is  above  the  threshold  range,  drying  is  relatively 
slow  for  the  SHM-derived  SWC  simulation.  In  this  case 
the  depth  of  the  shallow  surface  soil  layer  is  virtually 
irrelevant,  and  the  precise  value  of  SWC  is  of  little 
importance  except  insofar  as  it  was  above  the  threshold 
range.  Similarly,  when  Fr  is  large  (or  the  sky  is  overcast) 
the  threshold  range  or  substrate  drying  depth  is  irrele¬ 
vant.  For  that  reason,  the  increased  spatial  detail  in  the 
derived  Fr  may  yet  prove  to  be  more  useful  than  the 
SWC  values. 

These  results  underscore  the  importance  of  soil  type, 
which  is  a  controlling  factor  in  determining  the  hy¬ 
draulic  conductivity  and,  therefore,  the  threshold  drying 
range.  The  problem  in  predicting  decoupling,  however, 
may  remain  indeterminate,  because  the  composition  of 
soils  is  actually  a  mixture  of  different  types.  SWC  mea¬ 
surements  taken  at  five  locations  within  20  m  of  each 
other  at  an  individual  ARM  CART  station  can  exhibit 
variations  of  SWC  on  the  order  of  50%  of  each  other 
due  to  differential  drying  and  soil  types  across  the  site. 
It  is  tempting  to  believe  that  higher  spatial  resolution 
in  mesoscale  models  may  better  resolve  spatially  het¬ 
erogeneous  soil  layers,  but  the  variable  nature  of  soils 
and  soil  water  content  (on  the  order  of  1  m)  may  make 
this  goal  impossible.  Moreover,  it  is  not  obvious  how 
one  averages  soil  properties  in  such  an  inhomogeneous 
mixture,  and  a  similar  problem  also  presents  itself  in 
averaging  Fr.  At  the  heart  of  the  averaging  problem  is 
the  hugely  nonlinear  relationship  between  soil  hydraulic 
conductivity  and  soil  water  content. 

7.  Conclusions 

Although  we  have  not  demonstrated  that  a  derived 
SWC  is  necessarily  superior  to  a  default  value,  the  paper 
does  underscore  the  importance  of  rapid  surface  drying 
(decoupling)  in  prescribing  the  surface  energy  budget. 
Introduction  of  a  derived  Fr,  however,  does  appear  to 
be  useful  in  predicting  rapid  soil  drying.  Rapid  drying 
of  near-surface  soil  was  simulated  by  MM5/SHEELS  in 
many  locations  during  this  case  study,  particularly 
where  vegetation  cover  was  less  than  0.6.  Given  strong 
insolation,  and  a  vegetation  cover  less  than  0.6,  decou¬ 
pling  occurs  where  the  initial  soil  water  content  in  the 
top  10  cm  lies  within  a  threshold  range  (typically  be¬ 
tween  about  one-quarter  and  one-half  of  soil  saturation) 
determined  by  soil  type.  Based  on  simulations  for  the 
ARM  CART  SGP  site  in  central  Kansas  and  Oklahoma, 
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these  SWC  threshold  ranges  were  determined  for  var¬ 
ious  soil  types. 

Rapid  drying  was  accompanied  by  an  increase  in  sen¬ 
sible  heat  flux,  surface  skin  temperatures,  and  daily 
maximum  air  temperatures.  Although  decoupling  was 
evidenced  in  all  four  simulations  at  one  location  or  an¬ 
other,  the  SHM-derived  soil  moisture  remained  mostly 
above  the  threshold  range,  in  contrast  to  the  lower  de¬ 
fault  SWC  from  climatological  values.  The  most  ac¬ 
curate  temperature  forecasts  were  simulated  using  the 
derived  Fr  and  default  soil  moisture  fields.  Largest  errors 
in  maximum  temperature  were  found  for  the  derived  Fr 
and  the  derived  SWC,  which  contributed  to  more  evapo- 
transpiration  and  less  surface  heating  than  for  simula¬ 
tions  with  the  default  values. 

The  key  points  in  this  paper  about  the  simulation  of 
rapid  surface  drying  in  the  ARM  CART  SGP  site  cov¬ 
ering  central  Kansas  and  Oklahoma  are  as  follows. 

1)  Fractional  vegetation  cover  derived  from  biweekly, 
remotely  sensed  NDVI  data  provides  some  improve¬ 
ment  over  default  estimates. 

2)  Soil  moisture  simulated  from  a  hydrological  model 
with  finely  resolved  vertical  resolution  shows  signs 
of  decoupling  of  the  soil  layers  under  appropriate 
atmospheric  conditions.  However,  these  values  of 
soil  moisture  are  particularly  sensitive  to  the  soil 
parameterization  and  soil  type. 

3)  Higher-vertical  resolution  soil  models,  such  as  that 
used  in  SHEELS,  have  the  ability  to  simulate  rapid 
drying  of  the  upper  soil  layers  model  at  each  point. 

4)  Rapi9  soil  drying  (decoupling)  occurs  within  the 
0-2-cnl  layer  when  the  average  SWC  over  the  top 
10  cm  falls  into  certain  threshold  ranges  that  are 
dependent  on  soil  type. ' 

5)  The  effect  of  near-surface  SWC  and  soil  type  is  much 
larger  than  that  of  Fr,  with  the  exception  of  regions 
where  Fr  exceeds  80%  and  therefore  vegetation  and 
root-zone  soil  moisture  are  dominant.  More  accurate 
specification  of  vegetation  fraction  improved  the 
forecasts  overall,  but  the  derived  SWC  produced 
large  errors  in  simulated  temperatures,  possibly  be¬ 
cause  of  an  incorrect  interpolation  from  SHM  to 
SHEELS. 

6)  Future  modifications  to  SHEELS  will  ensure  that 
more  realistic  SWC  will  be  included.  The  accuracy 
of  such  models,  however,  may  be  limited  by  the  in¬ 
herent  local-scale  variability  of  soil  type  and  SWC 
and  the  highly  nonlinear  behavior  of  soil  hydraulic 
conductivity  as  a  function  of  SWC. 

Points  1-4  suggest  that  much  improvement  in  land 
surface  parameterization  can  be  accomplished  through 
remote  estimation  of  large-scale  biophysical  properties. 
Above  all,  it  is  necessary  to  include  a  high-resolution 
soil  model  to  capture  decoupling.  Increased  data  reso¬ 
lution  and  availability  will  make  it  even  easier  in  the 
future  to  simulate  soil  drying  at  smaller  scales,  as  new 


relationships  between  remotely  sensed  indices  and  bio¬ 
physical  variables  can  be  developed. 

Points  5  and  6  introduce  a  note  of  caution,  that  is, 
that  the  increased  complexity  of  land  surface  schemes 
and  measurements  does  not  always  result  in  improved 
predictions.  Soil  moisture  and  soil-type  data  are  so  cru¬ 
cial  for  predicting  decoupling  and  so  spatially  variable 
that  the  problem  of  predicting  decoupling  may  be  in¬ 
determinate.  An  important  goal,  nevertheless,  is  to  im¬ 
prove  the  land  surface  schemes  along  with  obtaining 
better  estimates  of  soil  and  vegetation  data. 
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